Ionic mixtures in two dimensions: from crystals to chain and dipole gases 
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The ground state of a two-dimensional ionic mixture composed of oppositely charged spheres is 
determined as a function of the size asymmetry by using a penalty method. The cascade of stable 
structures includes square, triangular and rhombic crystals as well as a dipolar pair gas and a gas of 
one-dimensional crystalline chains. Thereby we confirm the square structure, found experimentally 
on charged granulates, and predict new phases detectable in experiments on granular and colloidal 
matter. 
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Predicting the crystalline structures from first prin- 
ciples is one of the key problems in condensed mat- 
ter physics, material science, chemistry, geophysics and 
polymer physics [ij. For three-dimensional ionic crys- 
tals [1, it is common textbook knowledge Q that 
there are three basic equimolar structures whose stabil- 
ity depends on the ratio of the ion radii. For increas- 
ing asymmetry in the radii, the corresponding stability 
sequence involves the cesium-chloride, sodium-chloride, 
and zincblende structures. As indicated by the names, 
these structures are realized in nature for the (ionic) salt 
crystals CsCl, NaCl and ZnS respectively, but they 
show also up for a number of other equal-valency salt 
crystals and for mesoscopic oppositely charged colloidal 
particles which are suspended in a strongly deionized sol- 
vent [1, 0, • The stability is rationalized by minimizing 
the Madelung potential energy per particle pair with re- 
spect to various periodic equimolar candidate lattices at 
zero temperature and zero pressure. 

In two spatial dimensions, a similar problem arises 
by exploring the ground-state structure of ionic crystal 
monolayers in a model of oppositely charged disks with 
different radii. In fact, three different realizations of two- 
dimensional ionic crystals are possible on microscopic, 
mesoscopic and macroscopic length scales. First of all, 
crystalline sheets of molecular salts can be deposited on 
smooth substrates by using e.g. electrochemical meth- 
ods 8|. Secondly, oppositely charged colloidal particles 
0, |i| can be confined to a single layer, e.g. by using 
laser-optical trags or suspending them at a flat fluid- 
fluid interface [10|. Finally, there are granulate systems 
of millimeter-sized metallic balls which are oppositely 
charged and self-organize on a macroscopic plate into 



crystalline arrays [11|, . It is important to understand 
the different crystalline sub-structures both from a funda- 
mental point of view and for applications. For example, 
a control of the composite structures of colloidal crystals 
leads to new photonic [3l or phononic 14 1 band-gap ma- 
terials, to molecular-sieves 15[, to micro- and nano-filters 
with desired porosity [l6l | and to nanowires composed of 
individual particle strings [13] ■ It allows to steer protein 
crystallization 



Despite its fundamental importance, the stability of 
two-dimensional ionic monolayers has not been addressed 
for asymmetrically sized ions. Previous theoretical stud- 
ies have mainly focused on ionic criticality [l^l at finite 
temperature and on equal ion size [20| . In this Letter, 
we predict the ground state structures for two differ- 
ent relevant set-ups of oppositely charged spheres: in 
the first, all centers of mass belong to a common plane 
which corresponds to a situation where both species are 
confined in an interface [lij . In the second set-up, all 
particle surfaces touch a common plane corresponding to 
spheres in contact with the same plan ar substrate where 



they are confined by e.g. gravity [ll|,ll2]- Using a novel 
penalty method, a wealth of different stable structures is 
found. These include periodic crystalline arrays, on the 
one hand, with square, triangular and rhombic unit cells 
and with both touching and non-touching large spheres. 
On the other hand, a dipolar gas of particle pairs and a 
gas of particle chains with three different internal chain 
structures are also stable for strong size asymmetries in 
the substrate set-up. All these structures are detectable 
in experiments on granular and colloidal matter and in 
adsorbed crystalline layers of molecular salts. 

The model system used in our study are equimolar 1 : 1 
mixtures of large ions (component A) with radius Ra 
carrying a positive charge qa and small ions (component 
B) with radius Rb carrying a negative charge qs {\cia\ = 
I (7b I = q) such that the whole system is electroneutral. 
The size asymmetry < Rb/Ra < 1 is denoted with a. 
These constitutive ions interact via a pairwise potential 
composed by a Coulombic and hard-core part 



iir>Ra 
if r < Ra 



Rp 



{a,f3 = Aoi B) 



(1) 

where r is the center-center distance between the ions 
a and /3. The stability of the crystalline structure (at 
zero temperature) is ensured by steric interactions of the 
hardcore form. Our objective is to determine the sta- 
ble structures at zero pressure and zero temperature by 
numerical minimization of the total pot ential energy. 
Traditional minimization schemes [2lL typically re- 
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quire a continuous pairwise potential form in order to 
localize the minimum. The discontinuous hard-core po- 
tential splits the parameter space into various distinct 
regions which hampers straight-forward numerical mini- 
mization. For potentials which involve a hard-core part 
we have developed here a new technique to overcome this 
difficulty that relies on the so-called penalty method [2^ 
which was hitherto applied to geometric packing prob- 
lems of hard bodies. The key idea is to relax the non- 
overlap condition (i.e., r > Ra + Rp) by introducing an 
auxiliary penetrable pair interaction: 



+ /i max(i?a +Rp~r,0), (2) 



where /i > is a penalty parameter which is larger than 
the gradient of the potential at contact. This param- 
eters gives a finite energy penalty to any overlapping 
configuration. If fi is finite but sufficiently large, the 
total potential energy for the auxiliary potential Vaf3{r) 
has exactly the same minimal configuration as that for 
Uafsir)- However, numerically the potential energy land- 
scape is now continuous such that standard minimization 
routines like the simplex algorithm [3] can be applied. 
Though the penalty technique itself is applicable to any 
spatial dimensionality and hard particle shape, we ex- 
ploit it here to predict the ground state for binary hard 
charged spheres. 

In detail, we consider a parallelogram as a primitive 
cell which contains ua v4-particles and ub i?-particles. 
We restrict ourselves to the case ua < 3,n_B < 3. This 
cell can be described geometrically by the two spanning 
lattice vectors a and b. The position of a particle i (of 
species A) and that of a particle j (of species B) in the 
parallelogram is specified by the vectors rf = {xf,yf) 
and = {xf,yf), respectively. Thereby, the new po- 
tential energy function that needs to be minimized at 
zero pressure and zero temperature as a function of the 
crystalline lattice parameters reads 

t/total = 2 E E E' - + 

a=A,Bi,j = l L 

+ EEE^-4B(rf-rf +L), (3) 

i=l = l L 

with L = fca + where k and / are integers. The sums 
over L in Eqn.[3]run over all lattice cells where the prime 
indicates that for L = the terms with i = j are to be 
omitted. In order to handle efficiently the long range na- 
ture of the Coulomb interaction, we employed a Lekner- 
summation 25| . 

We now consider two different set-ups, the "interfa- 
cial model" and the "substrate model" . In the inter- 
facial model which can be considered as a purely two- 
dimensional situation, the centres of all spheres are con- 
fined to a plane, see Fig. [Ha) for a side view of a config- 
uration. In the substrate model, on the other hand, all 



spheres are touching the same underlying plane, see Fig. 
|3l[a) for a side view of a configuration. 

For the interfacial model, the stability phase diagram 
is shown versus the size asymmetry a in Fig. [ijb). By 
increasing cr, the following phase cascade occurs: 

A (touching, ^ 2) A (touching, iV^ = 3) 
A (non-touching, Nc ^ 3) ^ i?ft.(touching, iVc = 3) ^ 
□ (touching, Nc — 4:) ^ □(non-touching, Nc — 4) 

where the symbols A and □ stand for triangular and 
square unit cells of the big ions A, respectively, and Rh 
corresponds to rhombic unit cells, as illustrated by the 
top views of the crystalline structures in Fig. [Hb). The 
topological aspect of these four phases can be discussed in 
terms of contact between the large spheres. More specif- 
ically, a "touching" configuration involves connected big 
spheres while a "non-touching" one implies disconnected 
big spheres, see Fig. [UJb). In particular, the "touching" 
triangular and rhombic phases are both characterized by 
connected A-spheres in contrast to the "non-touching" 
A and □ structures which possess disconnected large 
spheres. Additionally, the ionic coordination number Nc 
defined by the number of A particles touching a single B 
particle, is another relevant topological characteristic. 

As far as the phase transitions reported in Fig. [T] 
are concerned, the following scenario takes place. At 
vanishing small ion size (cr — > 0), the small ions get 
squeezed between two big ions so that the three cen- 
ters of mass lie on a same line. The first transition 
A (touching, Nc = 2) ^ A (non-touching, Nc = 3) via 
the special structure A (touching, Nc = 3), character- 
ized by an increase of the number of contacts between A 
and B ions, occurs at cr = — 1. This special point 
corresponds to a compact triangular structure where a 
small ion has three contacts with neighbouring big ions, 
allowing a continuous transition. The second transition, 
A (non-touching, Nc = 3) ^ i?/i(touching, A'c = 3), oc- 
curring at cr = 0.297 is discontinuous as signaled by a 
jump of the angle between the two adjacent sides of the 
unit cell. The third transition i?/i(touching, A'^ = 3) ^ 
□ (non-touching, Nc = 4) via □(touching, Nc — 4), oc- 
curring at cr = V2 — 1 is continuous. A remarkable fea- 
ture, en passant, is the stability of the square phase over 
a wide range of the size ratio cr. 

For the two non-touching phases (A, □) one can assign 
two new Madelung constants M associated to the (lat- 
tice) electrostatic energy per molecule (i.e., a pair of ions 



A and B) E = -M „ . 
appropriate lattices provide: 



Lekner sums carried on the 



Ma = 1.542, Mo = 1.616 

for the non-touching triangular and square structures, 
respectively. As expected, this reported value of 
lies perfectly between that of the one dimensional lat- 
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FIG. 1; Stable structures of oppositely charged spheres versus 
their size asymmetry a — Rb/Ra in the interface model, 
where all sphere centers fall on the same plane: a) side-view, 
b) (scaled) energy per ion. The discontinuous transition is 
indicated by a solid bar. Continuous transitions are denoted 
by a broken bar. Unit cells of the corresponding stable phases 
are shown, where the big (small) have a radius Ra (Rb)- 



tice {Mid = 2 In 2) and that of a three-dimensional one 
(Mjvac; = 1.747). 

Next we focus on the area fraction covered by the 
spheres which is defined as (f) = '^^^a+^b) ^ ^fj{\}^ g^^^^ 
being the (projected) surface of the unit cell. The re- 
sults are sketched in Fig. [2l In the regime cr < — 1 
corresponding to compact triangular structures, there is 
enough space to host a small ion in the interstice offered 
by the touching big ions, so that the profile is identical 
to that of pure hard disk systems, where <j) — ^^^^ . In 
other words the location of the small sphere within the 
interstice does not alter the packing fraction. For non- 
touching triangular structures ("^ ^ 1 < < 0.297), 

d> varies like ^'^1+'^ ) _ jj^ f^i^^ rhombic phase regime 

(0.297 < a < \f2 — 1), the shape as well as the sur- 
face of the unit cell vary in a non-trivial manner. For 
large enough small ions, in the square phase regime, 
(\/2 - 1 < cr < 1), is given by 

We finally present results for the "substrate model" 
where all spheres are touching the same underlying sub- 
strate plane, see Fig. [3l^a) for a side-view. The stable 
crystalline structure and their energy per ion is shown in 
Fig. [2Ib) versus the diameter ratio a. For increasing size 
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FIG. 2: Area fraction as a function of the size aspect ratio 
cr = Rb/Ra- 

asymmetry a, a cascade of six structures is found: 

dipolar gas (non-touching, Nc = 1) ^ 
chain gas (non-touching, Nc = 1) ^ 
chain gas (touching, Nc~\)^ 
chain gas (touching, iVc = 2) ^ 
chain gas (non-touching, iVc = 2) ^ 
A (non-touching, Nc — 2>)^ □(non-touching, Nc = 4) 

In the limit of very large asymmetry (cr — > 0), dipoles 
perpendicular to the substrate plane are formed. Such 
parallel dipoles repel eachother and arrange into a crys- 
tal with diverging lattice constant and a coordination 
number Nc — 1. We call this state a "dipole gas" . When 
a increases, the dipole moments are getting gradually 
tilted relative to the substrate normal until a first-order 
transition towards a chain composed of dipolar pairs oc- 
curs where the individual dipoles are non-touching and 
the coordination numbers stays at Nc = 1. The dis- 
tance between neighbouring chains diverges since they 
are mutually repulsive. We are dealing therefore with a 
^^chain gas" at infinite dilution, i.e. the system is peri- 
odic in the direction along the string but with a diverg- 
ing interchain distance perpendicular to it bearing some 
analogy to smectic dipolar sheets [i^. Upon increasing 
a more, the chain structure changes continuously to an 
internal conformation with touching large spheres (still 
with Nc = 1). This structure then transforms continu- 
ously into a chain gas of non— touching large spheres with 
Nc = 2. Then, there are two crystalline structures ap- 
pearing known already from the interface model, namely 
a non-touching A and a non-touching □ lattice. The for- 
mer has only a tiny stability domain while the latter is 
stable along an enormous range of 0.517 < a < 1. We 
emphasize the striking emergence of the gas phases in the 
substrate model which are absent in the interface model 
and in three dimensions. 
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substrate model 



2Rg(side view) 
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FIG. 3: Stable structures of oppositely charged spheres ver- 
sus their size asymmetry g — Rb / Ra in the substrate model, 
where all sphere surfaces touch the same plane: a) side-view, 
b) (scaled) energy per ion. Discontinuous transitions between 
the structures are indicated by a solid bar. Continuous tran- 
sitions are denoted by a broken bar. Bottom views of the unit 
cells of the corresponding stable phases are displayed, where 
the big (small) have a radius Ra {Rb)- 



In conclusion, we have explored the stable ground state 
structures of two-dimensional ionic crystals for an "inter- 
face" and a "substrate" set-up at zero pressure. For a 1:1 
oppositely charged mixture of spheres with different di- 
ameters, various stable crystalline phases were identified 
including periodic crystals and chain and dipolar gases. 
Apart from adsorbed molecular salt systems, the struc- 
tures of the interface set-up can be verified in suspensions 
of oppositely charged particles. The substrate set-up, on 
the other hand, is realized for oppositely charged gran- 
ular matter on a plane. Recent experiments on oppo- 
sitely charged granular sheets with cr w 1 [ll|, [13] have 
indeed revealed a stable □(non-touching) configuration 
which is confirmed by our calculations. More experi- 
mental investigations on systems with higher size asym- 
metry are performable and could open the way to see 
the A (non-touching) and the predicted chain and dipo- 
lar gases. 

The penalty method can in principle be applied to any 
other potentials which involve hard-body parts, both in 
two and in three dimensions. It would be interesting 
to see the stability phase diagram for different mixtures 
as e.g. colloidal hard-spheres mixtures with or without 
added nonadsorbing polymers. 
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